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1.  INTRODUCTION. 


Mathematical  models  of  terminal  ballistics  processes  contain  uncertainties  that 
have  two  fundamentally  different  roots.  First,  one  has  the  usual  stochastic  errors  of 
observations  and  responses  of  materials.  These  errors  can  be  handled  by  conventional 
statistical  methods  based  on  estimated  probabilities,  although  in  ballistic  problems  the 
limits  of  application  for  these  methods  sometimes  must  be  stretched  due  to  missing  or 
scarce  data.  A  second  source  of  inaccuracy  is  the  lack  of  knowledge  about  processes  of 
threat/target  interaction.  Details  of  these  processes  are  difficult  to  observe,  they  are 
often  not  well  understood  and  not  represented  by  established  theories.  For  the 
treatment  of  this  type  of  inaccuracies  one  can  use  the  concepts  of  fuzzy  set  theory  that 
has  been  developed  to  quantify  judgement.  Fuzzy  sets  can  also  represent  inaccurate 
data,  in  particular  if  an  error  distribution  of  the  data  is  not  known  or  does  not  exist  (for 
instance,  in  intelligence  data  or  in  future  weapon  concepts).  In  the  present  report  we 
assume  such  a  representation,  that  is,  the  data  are  assumed  to  be  fuzzy  numbers  and  the 
mathematical  model  is  assumed  to  be  a  fuzzy  function.  We  then  formulate  a  fuzzy 
regression  problem  based  on  the  minimization  of  a  least  squares  objective  function. 

Special  types  of  fuzzy  regression  problems  have  been  treated  by  several  authors. 
Tanaka  et  al.,  [12]  and  [13],  developed  linear  programming  algorithms  for  regression 
with  fuzzy  linear  models  and  special  types  of  fuzzy  data.  The  algorithm  of  [12]  was 
applied  by  Heshmaty  and  Kandel  [9]  to  economic  forecasts.  Diamond  suggested  in  [6] 
and  [7]  a  least  squares  approach  and  in  [8]  a  maximum  likelihood  approach  to  fuzzy 
linear  regression  problems  with  special  fuzzy  data.  CelmipS  has  treated  regression  of 
nonlinear  models  with  a  slightly  more  general  type  of  data  in  the  special  cases  where 
either  only  the  data  [3],  or  only  the  model  [4]  is  fuzzy.  In  the  present  paper,  the  results 
of  [3]  and  [4]  are  extended  to  a  general  regression  problem,  where  the  model  equations 
are  nonlinear  and  the  model  as  well  as  the  data  are  fuzzy.  Some  restrictions  are 
imposed  on  the  generality  in  order  to  obtain  a  method  which  is  algorithmically  simple. 
The  restrictions  involve  a  parametrization  of  the  membership  functions  which  are 
assumed  to  be  conical  functions.  Under  this  restriction,  one  can  solve  the  model  fitting 
problem  with  the  aid  of  available  regression  software  for  implicitly  formulated 
constraints,  supplementing  the  software  only  with  normalization  procedures  for  some  of 
the  input. 

In  Section  2  we  develop  some  simple  concepts  and  properties  of  a  special  fuzzy 
point  space.  These  results  are  used  in  Section  3  to  derive  a  formalism  for  a  fuzzy 
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regression.  In  Section  4  the  fuzzy  regression  is  applied  to  the  task  of  developing  a  model 
for  the  bending  of  armor  plates  due  to  loads  by  pressures  from  explosions  of  bare 
explosives. 


2.  FUZZY  GEOMETRY. 

In  this  section,  we  provide  basic  definitions  of  fuzzy  set  theory  and  develop  simple 
concepts  of  analytic  geometry  in  a  fuzzy  environment.  The  purpose  of  the  section  is  to 
establish  practical  tools  for  model  fitting,  and  to  provide  a  geometrical  interpretation  of 
the  fitting  and  its  results.  Therefore,  the  exposition  is  restricted  to  a  special  type  of 
fuzzy  spaces  for  which  one  obtains  algorithms  that  are  easily  applicable  to  real  life 
problems. 

2.1.  Basic  Definitions. 

In  presenting  the  basic  definitions  we  follow  in  essence  Zimmermann  [16].  A 
classical  set  is  normally  defined  as  a  collection  of  elements  of  e£E.  Each  element  can 
either  belong  or  not  belong  to  a  set  HQE.  The  set  can  be  described  by  a  characteristic 
function  which  equals  unity  if  e  is  an  element  of  H  and  is  zero  if  e  does  not  belong  to  H. 
We  call  such  a  set  a  crisp  set.  A  generalization  of  the  concept  of  a  crisp  set  is  obtained 
by  allowing  the  characteristic  function  to  assume  values  that  are  not  restricted  to  zero 
and  unity.  To  distinguish  such  a  function  from  a  characteristic  function  it  is  called 
membership  function. 

Definition.  A  fuzzy  set  H  in  E  is  a  set  of  ordered  pairs 

H  =  {e  ,  pH{e)  |  eGE  }  .  (2.1) 

Pufe)  is  called  the  membership  function  or  grade  of  membership  of  e.  The  range  of 
Pn(e)  is  a  subset  of  nonnegative  numbers  whose  supremum  is  finite. 

In  this  paper,  we  assume  that  the  range  of  the  membership  function  is  the  interval 
[0,1]  and  that  elements  with  a  zero  degree  of  membership  are  not  included  in  H .  If  the 
membership  function  assumes  only  the  values  0  and  1,  then  the  set  H  is  non-fuzzy 
(crisp)  and  is  its  characteristic  function. 

The  concept  of  a  fuzzy  set  allows  one  to  express  uncertainty  in  a  manner  which  is 
different  from  probabilistic  concepts  and  interval  analysis.  Thus,  for  instance,  the 
concept  "about  5"  might  be  represented  by  a  fuzzy  number  5  whose  membership 
function  /i5(x)  has  a  positive  maximum  value  at  x=5  and  decreases  to  zero  as  |x—  5| 
increases.  Zimmermann  [16]  discusses  many  examples  of  such  concepts  and  their 
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representation  as  fuzzy  sets. 


2.2.  Fuzzy  Point  Space. 

A  fuzzy  point  in  Rn  is  a  fuzzy  set  of  points  XzRn.  Let  the  membership  value  of 
any  element  of  the  set  be  given  by  a  membership  function  p(X).  In  this  paper,  we  only 
consider  a  particular  type  of  fuzzy  points  which  are  characterized  by  fuzzy  position 
vectors  with  conical  membership  functions. 

Definition.  A  regular  (not  regular)  fuzzy  vector  A  in  Rn  is  the  pair 

X={A,Pa)  ,  (2.2) 

where  A€Rn.  and  PA  is  a  positive  definite  (semidefinite)  (nXn)-matrix.  A  is  called  the 
apex  of  A  and  PA  is  called  the  panderance  matrix  of  X. 

With  this  definition,  a  fuzzy  vector  in  an  n-dimensional  space  is  specified  by 
n  (n+3)/2  real  numbers.  Let  QA  be  the  Moore-Penrose  generalized  inverse  of  PA  and 
let  XA QRn  be  the  linear  space  spanned  by  the  columns  of  PA.  We  denote  by  ||  ||4  the 

following  elliptic  norm  of  the  distance  between  an  arbitrary  X€Rn  and  A: 

(  KX-zt)1,  •^•(X-A))1'*  .  X  X~  A  e  XA 

II*  All.,  - 1  ^  otherwise  .  (2  3) 

Definition.  The  conical  membership  function  of  A  is 

MaPO  =  1  ~  min  {  1.  ■  (2.4) 

The  fuzzy  point  A  is  the  set  of  all  vectors  X  with  the  property  f/A(X}>  0.  If  the 
panderance  matrix  PA  is  positive  definite,  then  the  boundary  of  the  support  of  pA  and 
the  level  surfaces  pA(X)= constant  are  hyperellipsoids  in  Rn  and  (2.4)  defines  a  cone. 
The  conical  membership  function  pA  is  normalized  to  0<HA<1,  and  a  fuzzy  point  with  a 
conical  membership  function  is  a  convex  fuzzy  set  [14].  If  the  panderance  matrix  PA  is 
semidefinite  then  the  cone  is  degenerate  in  Rn  and  is  a  regular  cone  in  the  subspace  A'^. 

In  Rv  the  conical  membership  function  is  the  triangular  function 

t*a{x)  =  1  -  min  {  1,  \x  -  a\  /  sa  }  ,  (2.5) 

where  sa='\/P^  is  the  spread  of  the  function.  In  Rn,  the  membership  function  of  each 
component  of  A  is  a  triangular  function  with  a  spread  that  equals  the  square  root  of  the 
corresponding  diagonal  element  of  PA.  The  triangular  function  is,  of  course,  the 
projection  of  the  cone  (2.4)  on  the  corresponding  axis.  Figure  1  gives  a  schematic 
illustration  of  the  conical  membership  function  of  a  fuzzy  point  (vector)  in  Rn. 


The  off-diagonal  elements  of  PA  indicate  possible  interactions  [15]  between  the 
components  of  A.  These  elements  can  be  conveniently  represented  by  dimensionless 
concordances  which  we  define  for  a  regular  fuzzy  point  as  follows 

cki  =  Pki  /  VPkkPa  .  k,i=l,...,n,  (2.6) 

where  and  are  elements  of  a  concordance  matrix  and  the  panderance  matrix  PA, 
respectively. 

We  introduce  a  structure  in  the  space  of  fuzzy  points  by  the  following  definition  of 
discord  and  grade  of  collocation  between  fuzzy  points. 

Definition.  The  discord  between  two  fuzzy  points  A  and  B  is 

D(X,B)~  min  max  {  ||jf  -  A\\a,  \\X  -  B\\b  }  .  (2.7) 

We  note  that  the  discord  is  not  a  distance  in  the  sense  of  functional  analysis, 
because  it  does  not  satisfy  the  triangle  inequality. 

If  the  panderance  matrix  of  one  of  the  vectors  in  eq.  (2.7),  say  PB,  approaches  a 
zero-matrix,  then  HB{X)  approaches  the  characteristic  function  of  the  apex  B.  That  is, 
the  corresponding  B  approaches  the  crisp  point  B.  The  corresponding  limit  of  the 
discord  is 


Jim  D(X,B)=\\B-A\\a  .  (2.8) 

B  — ♦  B 

Definition.  The  grade  of  collocation  between  two  fuzzy  points  A  and  B  is 

l(A,B)  =  max  min  {  pA(X),  HB(X)  }  .  (2.9) 

X€Rn 

The  grade  of  collocation  and  the  discord  are  related  by 

1{X,B)  =  1  -  min  {  1,  D{X,B)  }  .  (2.10) 


2.3.  Fuzzy  Functions. 

We  interpret  functions  as  mappings  of  a  domain  in  R„  to  a  range  in  Rr.  Usually, 
we  shall  assume  that  the  mapping  function  is  twice  continuously  differentiable  with 
respect  to  all  arguments  of  interest,  i.e.,  with  respect  to  XERn  and  with  respect  to  any 
free  parameters,  if  present.  If  some  of  the  arguments  are  fuzzy  numbers  or  fuzzy 
vectors,  then  the  image  of  a  crisp  or  fuzzy  point  in  R„  is  a  fuzzy  point  in  Rr.  For 
consistency  with  crisp  mappings,  we  shall  always  assume  that  the  image  of  the  apex  of 
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the  argument  is  obtained  by  replacing  in  the  function  all  fuzzy  arguments  by  their 
apexes.  That  is,  the  apex  of  a  fuzzy  point  is  mapped  onto  the  apex  of  the  image.  The 
panderance  matrix  of  the  image  can  be  computed  in  terms  of  the  panderance  matrices  of 
the  argument  and  parameters  by  formulas  which  are  derived  in  this  section. 


2.3.1.  Linear  Functions. 

We  first  consider  a  mapping  of  Rn  onto  itself  by  an  affine  coordinate 
transformation,  and  derive  a  panderance  propagation  formula  for  such  mappings  from 
the  requirement  that  the  membership  value  of  any  X6Rn  should  be  the  same  in  both 
coordinate  systems. 

Proposition  1.  Let  Z,  Z0,  X  and  A€Rn,  D  be  a  non-singular  (nXn)~ matrix,  and 

Z  =  L(X)=D  X  +  ZQ  .  (2.11) 


The  image  of  the  fuzzy  point  which  is  represented  by  the  fuzzy  vector 

X  -  {  A,  PA  } 


(2.12) 


is  the  fuzzy  point  (vector) 

B={B,Pb},  (2.13) 

with  the  apex 

B=D  A  +Z0  (2.14) 

and  the  panderance  matrix 

Pb=D  Pa  Dt  .  (2.15) 

Proof.  Equating  the  membership  values  of  an  XzRn  in  the  new  and  old 
coordinate  systems,  respectively,  one  obtains 


HB(  L(X) )  —  nA{X)  , 


or,  from  the  definition  (2.4), 


\\L(X)  -  B\\%  =  ||X-  A\\* 


(2.16) 


in  a  neighborhood  of  A.  If  X—A,  then  the  right  hand  side  of  (2.16)  vanishes  and  the  left 
hand  side  produces  eq.  (2.14).  Using  this  result,  we  express  the  left  hand  side  of  eq. 
(2.16)  by 
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(D  X  +  Z0  ~B)t  ■  Qb  •  (D-X  +  Z0 -B)  = 

~[D-(X-A)]t-Qb-[D-(X-A)\- 
=  (X-A)t  ■  DT  ■  QB  -  D  {X-A)  . 

Hence,  from  eq.  (2.16) 

dt  qb  d  =  qa  , 


(2.17) 

(2.18) 


or,  because  D  is  not  singular, 


Pn  =  D  P, 


q.e.d. 


The  extension  principle  for  fuzzy  functions  [15]  states  that,  if  Z  is  a  function  of 
L(X)  and  B—L{X),  then  one  should  have  the  relation 


wo-z-o^  ■  ,219) 

The  formulas  (2.14)  and  (2.15)  are  componentwise  consistent  with  this  principle  also  in 
cases  where  the  affine  mapping  (2.11)  is  not  a  coordinate  transformation.  The 
consistency  is  a  consequence  of  the  following  proposition: 

Proposition  2.  Let  Z  and  Z0eRr,  X  and  AeRn,  D  be  a  non-zero  (rXn)-matrix  of 
rank  r  and 

Z  =  L{X)=D  ■  X  +  Z0  .  (2.20) 

Then  the  image  B  of  a  regular  fuzzy  vector  JT  as  defined  by  eqs.  (2.14)  and  (2.15) 
satisfies  the  extension  principle  condition  (2.19). 

Proof.  We  formulate  the  right  hand  side  of  eq.  (2.19)  as  the  following  constrained 
minimization  problem: 

Minimize 

W  =  (X-A)t  PX'  (X-A) 

subject  to 

D  ■  (X-A)  =Z-B  . 

We  solve  (2.21)  using  a  Lagrange  multiplier  vector  k  and  defining  a  modified  objective 
function 

ft  =  \  {X  -  A)t  ■  P2l  ■  {X  -  A)  -  kT  ■  (  D  ■  {X  -  A)  -  Z  +  B)  .  (2.22) 

The  normal  equations  for  the  minimization  of  ft  are 
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PA'  {X-A)-DT  ■  k  =  0  , 

D  - {X- A) -Z  +  B  =  0.  ^2-23) 

Solving  the  first  eq.  (2.23)  for  k,  one  obtains 

k  =  (D  ■  PA  ■  DT)~l  ■  D  -  {X-A)  .  (2.24) 

We  substitute  this  expression  into  the  first  equation  (2.23)  and  multiply  the  result  from 
left  by  (X— A)t,  obtaining 

{X-A)t  PX1(X-A)=(X-A)t  Dt  (D  Pa  Dt)-1D  {X-A)  .  (2.25) 

This  equation  holds  at  the  solution  of  the  minimization  problem.  The  left  hand  side  of 
the  equation  is  the  sought  for  minimum  value  of  V$\  The  right  hand  side  can  be 
simplified  by  using  the  second  eq.  (2.23)  and  eq.  (2.15)  to  yield  (Z —B)T  ■  Pg1  (Z  —R). 
Hence, 

( (X  -  A)t  P?  ■  (X  -  A)  U  =  (Z  -  Bf  ■  P;'  ■  (Z  -  B)  ,  (2.26) 

q.  e.  d. 

Corollary  1.  If  the  mapping 

Z  =  L(X)  =  D  X  +  Z0  (2.27) 

maps  X  into  a  subspace. of  Rn  and  fiB(Z)  is  defined  by  eqs.  (2.14)  and  (2.15),  then  (2.19) 
holds  in  that  subspace. 


Proof.  Let  the  dimension  of  the  subspace  be  s.  Then,  by  a  proper  coordinate 
transformation,  (2.27)  can  be  rearranged  in  the  form 


Zy  —  ■  X  +  Z io  , 

Z2  =  0  , 


(2.28) 


where  D x  is  a  (sXn) -matrix  with  the  rank  s.  The  Proposition  2  applies  to  the  first  s 
equations  (2.28),  proving  the  corollary. 

If  A  is  not  a  regular  fuzzy  point  in  Rn  then  only  the  subspace  XA  (defined  in  the 
paragraph  following  eq.  (2.2.))  where  X  has  a  positive  membership  value  is  relevant  for 
the  mapping.  This  subspace  is  mapped  by  (2.20)  into  another  linear  subspace  ZACRn. 
Proposition  2  and  Corollary  1  directly  apply  to  the  mapping  of  XA  into  ZA. 


So  far  we  have  considered  a  crisp  linear  function,  eq.  (2.20),  with  a  possibly  fuzzy 
argument  X.  Now  we  let  the  function  be  fuzzy  and  assume  for  simplicity  that  the 
argument  X  is  crisp.  That  is,  we  consider  the  function 
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z  =  L(X)  =  D  X  +  Z0  ,  (2.29) 

where  some  of  the  elements  of  D  and  Z0  are  fuzzy  numbers  or  vectors.  By  definition,  a 
fuzzy  linear  function  is  a  subset  in  a  function  space  S.  That  space  consists  of  all  linear 
functions  which  map  Rn  into  Rr.  Its  dimension  is  r(n+l),  i.e.,  the  number  of  elements 
of  D  and  Z0.  Each  element  of  S  is  defined  by  a  particular  set  of  the  elements  of  D  and 
Zo,  which  in  turn  has  a  well  defined  non-negative  membership  value,  given  by  the  apexes 
and  panderances  of  D  and  Z0.  We  assign  to  the  corresponding  function  L(X)  the  same 
membership  value  thus  defining  the  fuzzy  function  set  E(X).  The  image  of  a  crisp  X  by 
eq.  (2.29)  is  fuzzy  due  to  the  fuzziness  of  the  function:  a  given  X  produces  several  images 
in  Rr,  each  with  a  corresponding  membership  value.  Also,  a  given  ZtRr  possibly  can  be 
the  image  of  several  points  XeRn  whereby  each  of  the  corresponding  mappings  has  a 
specific  membership  value.  According  to  the  extension  principle  for  fuzzy  functions  we 
assign  to  the  image  Z=^L(X)  the  supremum  of  these  values  as  its  membership  value. 

To  calculate  the  apex  B  and  the  panderance  matrix  of  the  image  B  ~t[A)  we  can 
invoke  the  Proposition  2,  because  eq  (2.29)  can  be  considered  as  a  crisp  linear  function 
(defined  by  the  crisp  elements  of  X)  of  the  fuzzy  arguments  D  and  Z0.  Let  Tx  and 
T2=^{T2,PT2}  be,  respectively,  the  vector  of  all  crisp  and  all  fuzzy  elements  in  D  and 
Z 0,  and  let  eq.  (2.29)  be  expressed  by 

Z  =  F(X,  I,,  f2)..  (2.30) 


Then, 

B  =F(A,  TV  f2) 

has  a  panderance  matrix  which  formally  is  given  by 

dF 


P  =M- 
B  dT , 


T  2 


dTo 


(2.31) 


(2.32) 


We  present  two  examples.  First,  let  z  be  a  scalar  function 


z  =  f(X ;  T)  =  g(X ;  C,z0)  =  XT  ■  C  +  z0 


and  let  T  be  the  fuzzy  vector 


II 

IE-h 

\c 

r°, 

,Pt 

Then,  for  X—A,  the  apex  of  the  function  value  is 


(2.33) 


(2.34) 
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=  6  =At  ■  C  + 


zo 


(2.35) 


and  its  panderance  matrix  is 


Pt={AT,l)  PT-  (?)  • 


(2.36) 


As  a  second  example,  assume  that  in  eq.  (2.29)  the  panderance  matrix  PD  of  the 
fuzzy  elements  of  D  is  diagonal,  and  that  all  concordances  between  elements  of  D  and 
Z0  are  zero.  The  formula  (2.32),  evaluated  at  X=A,  then  is 


rB  = 


EJW?  o 

t— 1 


J^P, Dii^i 

t— 1 


0 

0 


.EJW? 

l—l 


(2.37) 


2.3.2.  Nonlinear  Functions. 

We  consider  twice  differentiable  functions 


Z  =  F(W)  ,  (2.38) 

where  Z€Rr  and  WeRn.  Let  the  argument  W  be  a  fuzzy  vector  C .  Then  the  value  of 
Z  is  a  fuzzy  vector  B  with  the  apex 

B  =  F(C)  (2.39) 


and  a  panderance  matrix  which  may  be  approximately  computed  by  expanding  the 
function  F{W)  at  W=*C.  The  linear  part  of  the  expansion  is 

L(W)  =  F(C)  + -ITU  ■  (W  ~  C)  -  (2-4°) 


where  dF fdW  is  evaluated  at  W=C.  Applying  the  panderance  propagation  formula 
(2.15)  to  the  linearized  function  (2.40),  one  obtains  the  approximation 

(  \T 


Pb 


dF 

dW 


dF 

dW 


(2.41) 


For  the  intended  application  (mathematical  model  fitting)  it  is  convenient  to  split 
the  argument  W  of  a  function  into  two  sets:  observables  X  and  parameters  T,  that  is, 
to  express  model  functions  F(  W)  in  the  form 
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Z  =  F(X,T )  .  (2.42) 

Either  X  or  T,  or  both,  can  be  fuzzy  and  the  approximate  formula  (2.41)  may  be 

'*** 

applied  accordingly.  In  particular,  if  X=A  and  T=C,  and  the  panderance  matrix  of 
the  vectors  A  and  C  is 


Pw  = 


Pa 

P  CA 


r  AC 
Pr 


then  the  linearized  form  (2.41)  of  the  panderance  propagation  formula  is 


dF  p 

dF 

T,  dF  p 

dF' 

T,  dF  p 

dF) 

ax  Pa 

dX 

+  dT  Pc 

dT 

+  dXF*c 

dT  J 

.  M..p  . 

+  dT  ^CA 


dF 


dX 


(2.43) 


(2.44) 


We  consider  as  an  example  the  scalar  function 

z  =  f(X,T,z0)=XT  ■  T  +  z0  (2.45) 

with  the  fuzzy  arguments  X=K,  T=C  and  z0=t;.  The  function  can  be  expanded  at  the 
apex  (A,C,v)  of  its  argument  to  yield 

f(X,T,z0)  =  ATC  +  v  +  CT{X-A)+AT(T-C)+(z0-v)  + 

+  (t-c)t  (x-a)  .  (2  46) 

Using  the  approximate  formula  (2.44)  we  obtain 

PB**CT  Pa  C  +  AtPc  A  +  Pv+2AtPca‘C  +  2PvAC  +  2AtPc„  (2.47) 

The  approximation  comes  about  by  the  neglect  of  the  last  summand  in  eq.  (2.46),  i.e., 
the  nonlinear  part  of  the  function.  Its  contribution  generally  is  small  if  T  is  close  to  C 
and  X  is  close  to  A . 


The  advantages  of  the  approximate  calculation  of  the  panderance  of  a  function  by 
eq.  (2.41)  are  twofold.  First,  the  formula  is  simple  and  can  be  easily  included  in 
algorithms  using  fuzzy  arithmetic.  Second,  it  produces  function  values  which  have 
conical  membership  functions.  That  is,  with  this  approximation,  one  stays  in  conical 
fuzzy  vector  spaces  (including  degenerate  vectors  with  semidefinite  panderance 
matrices).  The  disadvantage  is,  as  in  any  approximation,  that  in  certain  cases  the 
approximation  may  not  be  sufficiently  accurate. 
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2  3.3.  Fuzzy  Function  Space. 

In  analogy  to  the  linear  functions  discussed  above  we  may  consider  the  nonlinear 
function 

Z  =  F{X,T)  (2.48) 

with  a  fuzzy  parameter  vector  T  as  a  fuzzy  function  of  a  crisp  argument: 

Z  =  F(X)  .  (2.49) 

To  define  the  fuzzy  set  of  functions,  F ,  we  specify  a  corresponding  membership  function. 
In  the  case  F(X)=F(X,T),  each  member  of  the  set  F  is  given  by  a  particular  parameter 
value  te  T.  We  assign  to  the  function  F(X,t)  the  membership  value  of  t: 

.  (2.50) 

In  this  paper,  we  only  consider  fuzzy  functions  of  the  described  type,  that  is,  we  assume 
that  a  fuzzy  function  always  is  defined  in  terms  of  a  fuzzy  function  parameter. 

To  measure  the  separation  between  elements  of  a  fuzzy  function  space,  we  define 
the  following  norm  ||  •  ||fr  for  the  distance  of  the  crisp  function  F(X,t)  from  the  crisp 
"apex"  function  F{X,  T)  with  the  membership  value  one: 

II F(X,t)  -  F(X,T)\\ft  =  lit  -  T|lr  -  [  (<  -  T)T  ■  Pjl  •  (t  -  T)  ]1/2 .  (2.51) 

(If  T  is  not  a  regular  point,  then  Pf1  is  replaced  by  QT  and  the  norm  defined  in 
accordance  with  eq.  (2.3)  ).  With  this  definition,  one  has  in  analogy  to  (2.4) 

nM  F(X,t) )  =  1  -  min  {  1,  ||F(X,f)  -  F(X,T)  ||fT  }  .  (2.52) 

The  discord  and  the  grade  of  collocation  between  two  fuzzy  elements  of  the  fuzzy 
function  space  are  defined  in  analogy  to  eqs.  (2.7)  and  (2.9),  respectively. 

2.4.  F uzzy  Equations. 

A  fuzzy  equation  is  a  fuzzy  set  of  equations.  One  obtains  a  fuzzy  equation,  for 
instance,  by  setting  a  fuzzy  function  equal  to  a  crisp  or  fuzzy  zero.  Let  X6Rn,  T€Rp 
and  the  function  F(X,T)ERr.  We  further  assume  that  F  is  at  least  twice  continuously 
differentiable  in  a  proper  domain, 

r  <  n  and  r  <  p  .  (2.53) 

We  define  a  fuzzy  function  as  in  Section  2.3.3.,  and  obtain  a  fuzzy  equation  by  setting 
the  function  equal  to  a  crisp  vector  with  r  vanishing  components  : 
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F(X)  =  F{X,  T )  -  0  .  (2.54) 

We  shall  assume  that  the  right  hand  side  of  a  fuzzy  equation  always  is  a  crisp  vector. 
By  permitting  it  to  be  a  fuzzy  vector  one  does  not  gain  additional  generality,  but  merely 
increases  the  number  of  fuzzy  parameters  by  r.  That  is,  the  fuzzy  equations 

F(X)  =0  (2.55) 

and 

G(X)=F(X)  -0=0  (2.56) 

are  equivalent  sets  of  equations  in  a  sense  which  we  shall  describe  shortly.  Any  crisp 

■**w  _ 

T  defines  a  particular  crisp  equation  of  the  fuzzy  set  (2.54).  That  crisp  equation 
either  has  no  solution  or  defines  a  crisp  solution  Xp(t).  We  assign  to  the  equation  and  to 
its  solution  the  membership  value  of  the  function  F(X,t),  that  is,  the  value  nT{t). 

The  solution  Xp(t)  of  the  crisp  equation  F(A,i)=0  is  in  Rn  a  (n  —  r)-dimensional 

hypersurface  which  is  the  set  (locus)  of  all  points  X  that  satisfy  the  equation.  The  fuzzy 

equation  (2.54)  defines  in  Rn  a  fuzzy  surface  X  F,  that  is,  a  fuzzy  set  of  (n  —  r)- 

dimensional  surfaces.  The  surface  corresponding  to  the  parameter  value  t  is  assigned 

the  membership  value  /iT(f).  The  fuzzy  surface  XF  contains  all  points  X  which  satisfy 

*>** 

eq.  (2.54)  for  some  f€7\  To  an  X€XF  we  assign  the  largest  membership  value  of  all 
solution  surfaces  that  pass  through  X,  that  is, 

HxAX)  -  t  Fs{$>t)_0  Pni  HX,t) )  =  t  Fsupf)_0  Hj(t)  •  (2.57) 

If  X$XF,  i.e. ,  if  X  does  not  solve  the  eq.  (2.54)  for  any  t£  T,  then  we  assign  to  X  the 
membership  value  zero. 

Now  we  consider  eqs.  (2.55)  and  (2.56).  Their  equivalency  is  related  to  the 
following  proposition. 

Proposition  S.  Let  X6 Rn,  C€Rrl  C  ={C,Pc},  and  the  function  G(X)€Rr.  Then 
the  solutions  of 

G(X)  =  C  (2.58) 

and 

F{X,C)  =  G{X)-C  =  0  (2.59) 

are  identical  fuzzy  sets. 

Proof.  Eq.  (2.58)  is  a  fuzzy  set  of  equations.  Let  its  solution  be  X G.  Because  one 
side  of  the  equation  is  crisp,  one  may  use  eq.  (2.57)  to  construct  the  solution  in  terms  of 
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its  membership  value.  The  result  is 


Vxg(X)  =  t  su 

t  :  Cr  I 


(2.60) 


On  the  other  hand,  according  to  eq.  (2.57)  the  solution  of  eq.  (2.59)  is  a  fuzzy  set  with 
the  membership  function 


MxA*)  -  t  MO  •  (2-61) 

Because  the  delimiters  for  the  maximization  of  /i^l)  in  eqs.  (2.60)  and  (2.61)  are 
equivalent  crisp  expressions,  pXG(X)=nXF(X),  q.e.d. 

Corollary  2.  The  solutions  of  eqs.  (2.55)  and  (2.56)  are  identical. 

Proof.  This  is  the  special  case  C  =  0  of  Proposition  S. 

Proposition  3  and  Corollary  2  show  that  certain  algebraic  manipulations  of  fuzzy 
equations  are  permitted  in  the  sense  that  they  do  not  affect  the  solution  of  the 
equations.  We  note,  however,  that  not  all  manipulations  which  are  permitted  for  crisp 
equations  can  be  applied  to  fuzzy  equations.  For  the  present  paper,  the  important 
consequence  of  Proposition  S  is  that  one  can  assume  on  the  right  hand  side  of  a  fuzzy 
equation  a  crisp  vector  with  zero-valued  components. 

The  definition  (2.57)  of  the  membership  function  of  XF  is  not  practical  for 
numerical  algorithms.  We  provide,  therefore,  another  equivalent  definition  in  terms  of  a 
norm  for  the  distance  between  X  and  the  solution  of  F(X,T)=  0,  i.e.,  between  any 
X£Rn  and  the  ctisp  solution  surface  XF(T).  We  define  such  a  norm  by 

\\X-XF\\XF=  inf  ||l  —  T\\t  (2.62) 

t :  F(X,t)—0 

and  assign  to  ||X— Xyr||  an  infinite  value  if  F(X,t)^0  for  all  teRp.  With  this  definition, 
eq.  (2.57)  is  equivalent  to 

f*xA^ 0  —  1  ~  min  {  1,  IIX  —  XF\\XF  }  •  (2.63) 

This  formulation  is  better  suited  for  numerical  calculations  than  (2.57),  because  the 
norm  (2.62)  can  be  approximately  computed  as  shown  by  the  following  proposition. 

Proposition  4.  Let  the  r  components  of  F(X,t)  be  linearly  independent  functions 
of  t.  Then  the  following  approximation  of  (2.62)  is  exact  for  functions  F(X,t)  that  are 
linear  with  respect  to  t: 


-13- 


Wx-X, 


F»XF 


F(x,t y 


dF 

dt 


dF 


dt 


-1 


F(X,t) 


1/2 


(2.64) 


Proof.  We  use  the  definition  (2.62)  and  formulate  the  left  hand  side  of  eq.  (2.64) 
as  the  following  constrained  minimization  problem: 


Minimize 


subject  to 


W  =  {t  —  T)t  ■  Pf1  •  (t  —  T) 


F(X,t)=  0  . 


(2.65) 


Using  a  Lagrange  multiplier  vector  k,  one  obtains  for  the  problem  (2.65)  a 
modified  objective  function 

ft  =  ±{t  -  T)t  ■  P f1  •  [t  -  T)  -  kT  ■  F(X,t)  .  (2.66) 

The  corresponding  normal  equations  are 


Pf 1  (t-T)-Fl  ■  k  =  0 


F(X,t)=  0  , 


(2.67) 


where  Ff—dF/dt.  After  elimination  of  k  from  the  first  equation  (2.67)  and 
multiplication  of  the  result  from  left  by  (t  —  T)t ,  one  obtains  the  following  formula  for 
||  t  —  r||2  at  the  solution: 

(t  -  T)t  ■  Pf 1  (t-T)=(t-  T)t  ■  Fj  ■  {F,  ■  PT  ■  FI)'1  ■  F,  (t-T)  .  (2.68) 
On  the  right  hand  side  of  eq.  (2.68)  we  approximate 

•  (t  -  T)  *  -F{X,  T)  (2.69) 

(because  F(X,t)—0  ),  and  also  use  the  argument  (X,T)  instead  of  (X,t)  in  the  product 
Ft-PT-Fl.  The  result  is  the  right  hand  side  of  eq.  (2.64). 

If  F  is  linear  with  respect  to  t,  then  (2.69)  exactly  holds,  and  Ft  is  independent  of 
t.  Therefore,  (2.64)  is  exact  in  the  linear  case,  q.e.d. 

Now  we  investigate  the  structure  of  the  membership  function  iiXF{X)  in  the 
vicinity  of  the  apex  of  the  function.  The  apex  itself  is  the  crisp  hypersurface  XF(T), 
defined  by  F{X,T)=  0.  The  level  surfaces  or  IIX— ^11^=1  — 7,  approximately 
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are  given  by  the  equation  (see  (2.64)  ) 

Ft  ■  (Ft  ■  PT  ■  Fj)~x  ■  F  =  (1  -  7)2  •  (2.70) 

Let  YeXF{  T),  i.e.,  Y  be  a  point  of  the  apex  of  Xp{T).  Then,  F(Y,T)=  0  and,  expanding 
F  in  eq.  (2.70)  in  terms  of  X  —  Y  and  only  keeping  the  first  non-zero  term,  one  obtains 
for  the  level  surface  the  equation 

(X-Y)T  F$  (Fi  PT  F?)-1Fx  (X-Y)=(l-if  ,  (2.71) 

where  Fx=dF  /dX.  Within  the  accuracy  provided  by  the  one-term  approximation,  the 
structure  of  the  level  surface  in  the  vicinity  of  Y  is  characterized  by  the  matrix 

QY  =  r$(F,PTFjr'Fx.  (2.72) 

We  recall  that  F€Rr  and  r<n.  The  rank  s  of  the  matrix  QY  is,  therefore,  at  most 
equal  to  r<n.  (It  equals  r  if  Fx  has  a  full  row  rank.)  If  s=n,  then  the  solution  X  F  is  a 
fuzzy  point  in  Rn,  and  the  surfaces  ||X— Xf||^jp=l— ■ 7  are  hyperellipsoids  in  Rn.  If 
s<n,  then  the  matrix  Qy  is  only  semidefinite  and  the  surfaces  are  hypercylinders.  If  7 
is  close  to  one,  then  the  boundaries  of  the  support  of  fiXF(X)—~i  are  in  case  s<n 
cylindrical  surfaces  that  are  approximately  parallel  to  the  apex  surface  XF(T). 

If  the  apex  XF(T)  is  a  point  in  Rn  then  its  distance  norm  ||  •  ||A  from  a  fuzzy  point 
A  is  defined  by  eqs.  (2.2)  to  (2.3).  If  the  apex  X^(  T)  is  a  surface  in  Rn  then  we  define 
the  distance  by 

\\XF  -  /111 A  =  min  [  (X  -  A)t  ■  P?  ■  [X  -  A)  ]'/2  = 

X  €Xp 

(2.73) 

=  min  [(AT-A)r  •  Pj1  •  (X-A)]1/2  . 

X :  F(X,  r)—0 

An  approximate  formula  for  llA'jp— A\\a  is  provided  by  the  following  proposition. 

Proposition  5.  Let  the  r  components  of  F  be  linearly  independent  functions  of  X. 
Then  the  following  approximate  formula  for  (2.73)  is  exact  for  functions  F(X,t)  that  are 
linear  with  respect  to  X: 

IUf  -  M\a  *  (  F T  ■  iFx  PA  ■  F]c)-'  F  }'n  ,  (2.74) 

where  F  and  dF /dX  are  evaluated  at  {A,T). 

Proof.  The  proof  follows  the  same  steps  as  the  proof  of  Proposition  4,  with 
appropriate  changes  of  the  arguments  of  F(X,  T). 

The  separation  between  a  fuzzy  point  A  and  the  fuzzy  surface  XF  may  be 
measured  in  analogy  to  (2.7)  by  the  discord 
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D(X,XF)  =  min  max  {  ||X  —  All^  ,  ||X  —  XF^\XF  }  (2-75) 

X  Ei?* 

The  grade  of  collocation  between  ^  and  X F  is  (see  eq.  (2.10) ) 

l(X,XF)  =  1  -  min  {  1  ,  D(X,XF)  }  .  (2.76) 

In  some  important  special  cases  the  discord  D(A,XF)  can  be  explicitly  computed. 
These  cases  are  covered  by  the  following  Propositions  6  and  7. 

Proposition  6.  Let  Fx  and  Ft  be  full  rank  constant  coefficient  matrices, 

F(X,t)  =  F0  +  FxX  +  Ft- t  (2.77) 

Aix^n  pA  be  non-singular, 

II*  -  A\\a  =  ( (X  -  A)t  ■  PAl  ■  (X  -  A)  I1/2  ,  (2.78) 

II*  -  XfWxF  =  I  nx,  T)T  •  (Ff  •  PT  •  •  F(X,  T)  }''2  ,  (2.79) 

and  L(X)  be  the  locus  of  points  where  \\X—A\\A=constant  is  tangent  to 
||  X— XF^\XF=*  constant. 

Necessary  and  sufficient  for  the  locus  to  be  a  straight  line  is  that 

f)  F,  PT  Ff=Fx  PA  Frx  (2.80) 

for  some  positive  0. 

Proof.  We  first  note  that  L(X)  by  definition  contains  the  point  A,  and  that  along 
any  straight  line  through  A  the  planes  tangent  to  \\X—A\\a— constant  are  parallel  to 
each  other.  The  tangent  plane  which  passes  through  a  point  C^A  is  spanned  by  vectors 
X  —  C  that  satisfy  the  equation 

(C  —  A)T  ■  PAl  ■  (X  —  C)  =0  .  (2.81) 

The  plane  tangent  to  \\X—XF\\XF= constant  which  passes  through  a  point  C$  XF(t)  is 
spanned  by  vectors  X—C  that  satisfy  the  equation 

F(C,T)t  ■  (Ft  ■  PT  ■  Ff)~l  ■  Fx  (X-C)=  0  .  (2.82) 

Let  B€Xf{T),  i.e. ,  F(B,T)—  0.  Then,  F(C,T)=FX  (C—B)  and  (2.82)  can  be  expressed 
in  the  form 

( C-Bf  ■  F%  ■  (F,  ■  PT  •  F,r)->  •  Fx  (X-C)  =  0  .  (2.83) 

Eq.  (2.83)  shows  that  planes  through  C  and  tangent  to  \\X—XF\\XF= constant  are 
parallel  to  each  other  for  all  such  C  which  are  located  on  any  straight  line  that  passes 
through  a  BtXF,  i.e.,  through  the  apex  of  X F. 
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We  find  a  point  BeXf(T)  which  also  is  a  point  of  the  locus  L(X)  by  solving  the 
following  minimization  problem: 

Minimize 

W  =  {B  —  A)t  ■  Pj1  ■  ( B-A ) 

subject  to 

F{B,T)=  0  . 

Using  a  Lagrange  multiplier  vector  k,  one  obtains  a  modified  objective  function 
W  =  ±  ( B-A)t  ■  PJ1  •  (B-A)  -  kT  •  F(B,  T) 
and  the  corresponding  normal  equations 

Pxl  (B-A)-F^k=  0 
F{B,T)  =  0  . 


(2.84) 


Eliminating  k  from  the  first  equation  (2.84)  one  obtains 

B-A=PaFI(Fx  Pa-  f £)->  '  Fx  ■  (B-A)  , 

or,  because  F(B,  T )  =  F(A,T)  +  Fx  ■  ( B—A )  =  0, 

B  —  A  =  —PA  Fx  (Fx  ■  PA  ■  FrxY'  ■  F(A ,  T )  (2.85) 

The  plane  tangent  to  \\X—A\\a=  constant  through  the  point  B  is  spanned  by  vectors 
X—B  which  satisfy  the  equation  (see  (2.81)  ) 

(B-A)T  PA'  (X-B)-O  , 

or,  using  eq.  (2.85), 

F(A,T)t(FxPaFZ)-'Fx(X-B)=  0  (2.86) 

On  the  other  hand,  the  surface  ||X— XF\\ XF= constant  which  passes  through  the  point  ,4 
has  a  tangent  plane  that  is  spanned  by  vectors  X— A  satisfying  the  equation  (see  (2.82)  ) 

F(A,T)t  ■  (Ft  ■  PT  ■  Fir1  ■  Fx  ■  (X  -  A)  =  0  .  (2.87) 

Necessary  and  sufficient  for  L(X)  to  be  a  straight  line  is  that  the  planes  defined  by  eqs. 
(2.86)  and  (2.87),  respectively,  are  parallel  to  each  other.  Inspecting  the  equations,  one 
sees  that  necessary  and  sufficient  for  the  planes  to  be  parallel  is  that  (2.80)  holds,  q.e.d. 

Proposition  7.  Let  the  premises  of  Proposition  6  hold,  and  (2.80)  be  satisfied. 

Then, 
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(2.88) 


nts"  v  \  -  a||x  •  IU  -  xF\\XF 

D[A'Xf)~  \\Xf-A\\a  +  \\A-Xf\\xf  ’ 

where  ||Xf— ^4(1^  is  given  by  eq.  (2.74)  and  ||/i— Xf|lxF  is  given  by  eq.  (2.79). 

Proof.  The  discord  D(A,X F)  is  defined  by  eq.  (2.75).  Because  || — ^4 1| A  as  well 
as  \\X—Xf\\xf  are  unbounded  continuous  functions,  differentiable  everywhere,  except  at 
their  zeros,  the  minimum  in  eq.  (2.75)  is  assumed  at  a  point  C  where 
II C—  A||a=||  C—Xf\\xf,  and  the  planes  tangent  to  \[X—A\[a= constant  and  tangent  to 
||X— Xjr|| xF=constant  coincide.  C  necessarily  is  a  point  of  the  locus  L(X),  or, 
according  to  Proposition  6  a  point  of  the  straight  line  through  A  and  B,  the  latter  point 
being  given  by  eq.  (2.85).  Therefore,  one  can  compute  D  by  first  finding  that  point  C 
between  A  and  B  where  both  norms  are  equal,  and  then  computing  the  value  of  the 
norms  at  X=C. 

Let 

C  =  A+<x(B-A)  (2.89) 

be  a  point  of  L( X).  Then,  by  definition 

\\C-A\\A=\a\  \\B-A\\A=\a\  .\\Xf-A\\a 

and 

II C  —  Xf\\xf  =  1 1  —  cv|  •  \\A  —  Xf\\xf  ■ 

Both  expressions  are  equal  for 

_ ll-A  —  Xf\\xf _ 

l|X/p  ~  -^ll/i  +  11^  —  -X/rll  xf 

Hence, 

D{X,XF)  =  ||  C  -  A\\a  =  a  •  ||Xj,  -  A\\a 

is  given  by  eq.  (2.88),  q.e.d. 

Corollary  S.  Let  F(X,t)=F0  +  FxX  +  F,  t  be  a  scalar  linear  function  of  X  and 
t.  Then  the  discord  D(A,XF)  is  given  by  eq.  (2.88). 

Proof.  If  F  is  a  scalar  linear  function,  then  eq.  (2.80)  is  satisfied  and 
Proposition  7  holds. 

Corollary  4 ■  Let  A  and  B  be  fuzzy  points  in  Rn,  and  P^^Pb-  Then  the 
discord  between  A  and  B  is  given  by 
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(2.90) 


Proof.  We  define  B  as  the  solution  of  the  fuzzy  equation 

F{X,B)  =  X-B  =0  . 


Propositions  6  and  7  hold  for  this  function  and,  therefore 


D{X,B) 


0\\B~A\\2a 

(1+0)\\B-A\\a 


0 

1+0 


■\\b-a\\a 


q.e.d. 


If  F(X,T)  is  not  linear  in  X  and  T,  then  the  discord  formula  (2.88)  may  be  applied 
to  a  linearized  approximation  of  the  equation  F[X,T)=  0,  e.g.,  to  a  fuzzy  set  of 
equations  with  the  elements 


F(X,t)  =  Fx  •  (X  -  V)  +  ,Ft  ■  (t  -  T)  =0  , 

where  Y€Xp{T).  For  scalar  F  this  yields  a  reasonable  value  for  the  discord  if  t  is 
sufficiently  close  to  T  and  X  is  sufficiently  close  to  XF(T).  If  F  is  not  scalar,  then  also 
the  condition  (2.80)  needs  to  be  satisfied  to  have  a  reasonable  approximation  of  the 
discord.  In  particular,  in  the  case  r=n  the  condition  (2.80)  is  important. 

We  make  a  final  remark  about  the  fuzziness  of  the  surface  F(X,T)=  0.  Its  spread 
is  defined  by  eqs.  (2.62)  and  (2.63).  It  is  easy  to  see  that  the  spread  depends  on  the 
formulation  of  the  equation  F(X,T)—0.  Therefore,  two  forms  of  the  equation  which 
produce  the  same  crisp  apex  surface  Xp(T)  (defined  by  F(X,T)=  0)  in  general  have 
different  spreads,  i.e.,  they  define  different  fuzzy  surfaces  X  F.  This  dependence  of  X F  on 
the  equation  formulation  is  an  intrinsic  property  of  fuzzy  equations. 


3.  FUZZY  MODEL  FITTING. 


3.1.  Fuzzy  Models. 

We  consider  mathematical  models  of  observable  events  that  are  formulated  as 
twice  differentiable  fuzzy  equations  among  components  of  a  fuzzy  observable  vector  X . 
That  is,  the  models  are  fuzzy  equations  of  the  type 

F{X,¥)  =0  ,  (3.1) 

where  F  is  a  twice  differentiable  vector  function.  It  is  convenient  to  assume  that  the 
model  parameter  vector  T  consists  of  two  parts,  T j  and  Tp.  The  first  part,  T j, 
represents  free  model  parameters.  The  second  part,  T  p,  consists  of  prescribed 
parameters.  (Either  of  these  can  be  an  empty  set,  of  course).  The  apex  T j  of  T  j  is 
determined  by  fitting  the  model  to  observations  while  the  apex  Tp  is  predetermined. 
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The  corresponding  panderance  matrices,  PTj  and  PTp,  of  these  vectors  may  be 
prescribed  or  determined  concurrently  with  other  unknowns  of  the  regression  problem. 

We  assume  throughout  the  rest  of  this  section  that  the  dimensions  of  F,  X  and  T 

are  r,  n  and  p,  respectively.  We  will,  however,  permit  the  functions  F{  and  the 

corresponding  dimensions  r{  and  n{  to  be  different  for  different  observation  vectors  Xt. 
For  simplicity,  we  assume  that  the  model  parameter  dimension  p  is  the  same  for  all 
observations.  This  is  not  a  restriction,  because  not  all  model  functions  need  to 
depend  explicitly  on  all  parameters. 

3.2.  Model  Fitting  Principles. 

We  determine  the  unknown  model  parameter  vector  T  by  maximizing  the 
membership  values  of  the  corrected  observations  as  well  as  the  membership  values  of 

the  model  surface  at  the  corrected  observations.  Because  the  membership  functions 

have  a  normalized  maximum  of  one,  the  maximization  is  equivalent  to  a  minimization 
of  the  deviations  of  the  membership  values  from  one.  In  particular,  we  minimize  the 
sum  of  the  squares  of  such  deviations. 

The  constraints  of  the  model  fitting  problem  are  the  equations 

F i(Xi  +  cit  T  +  e,)  =  0,  *  =  1 . s  .  (3.2) 

That  is,  local  adjustments  are  applied  to  the  apexes  X<  of  the  observations  Xt={X,,Px,} 
as  well  as  to  the  apex  Tt=(TJ,TJ)  of  the  model  parameter  vector.  The  local 
corrections  are  only  applied  to  the  fuzzy  components  of  X,-  and  T.  That  is,  if  a 
component  of  these  vectors  is  crisp,  then  the  corresponding  local  adjustment  is  zero. 
The  non-zero  adjustments  are  obtained  by  minimizing  the  objective  function 

a 

»'=E|(i-(<y,«  +  cj))2+(i-Mr  +  «,'))!]  (3.3) 

t=l 

subject  ot  the  constraints  (3.2). 

The  solution  of  the  model  fitting  problem  consists  of  the  model  parameter  vector 
T  and  the  set  of  the  local  adjustments  (residuals)  ct  and  e,-.  In  the  r?,-dimensional  space 
of  the  *-th  observation  set,  the  fitted  surface  (the  model)  is  the  fuzzy  solution  of  the 
equation  Ft(X,  T)=0.  It  is  reasonable  to  require  that  the  observations  X,  are  close  to 
the  fitted  surface.  In  our  fuzzy  environment  we  express  such  a  requirement  in  terms  of 
the  grade  of  collocation  between  X,  and  XFi,  that  is,  by 

l(Xi,XFi)  >  7,  ,  i  =  (3.4) 

The  choice  of  7,  depends  on  the  particular  application  and  expresses  the  modeller’s 
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confidence  in  the  adequacy  of  the  model  and  in  the  reliability  of  the  data  panderance 
estimates.  Typically  one  specifies  different  desired  values  of  "y ,•  in  the  course  of  an 
investigation  or  sets  all  desired  values  equal  zero  and  calculates  a  posteriori  the  attained 
grades  of  collocation.  The  attained  grades  then  might  be  used  as  outlier  indicators. 


Using  the  relation  (2.76)  between  ^(X^X Fi)  and  the  discord  D(XitX Fi),  and  the 
discord  formula  (2.88),  this  condition  can  be  approximated  by 


II*r  -  *JU  ll*i  -  *nllw,  ..  • 

Il-X/V  —  -^illjfi  +  II-X»‘  —  ^FilljCFi 

(3.5) 

or,  in  case  of 

a  scalar  linearizable  model,  by 

Ikll*.  ' Ml  r  , 

Vlklljs+vIMr  " 

(3.6) 

where  the  norms  are  given  by  eq.  (2.3)  in  terms  of  the  panderances  PXi  and  PT, 

Fi(Xi,T  +  ei) 

■’  FIXt.T) 

(3.7) 

and 

.  FM  +  chT) 

nx„r,  ' 

(3.8) 

The  corrector  e,-  appears  in  eq.  (3.6)  because  II  c»ll av  *s  not  norm  of  the  whole 
distance  of  the  observed  point  Xi  from  the  apex  X*,  but  only  of  a  part  of  it.  A  linear 
approximation  yields  llc»lljtt— -^.ll*,-  and,  correspondingly, 
II  e,ll  r~&i'  ll-^i  “ -X/’.-ll  XFi 

If  at  the  solution  all  the  nXi  and  *n  t^ie  expression  (3.3)  are  positive  then  one 
can  use  the  definition  (2.4)  and  express  the  objective  function  (3.3)  in  terms  of  the  norms 
II  cillx«  anc^  Ikll  T  as  follows: 

8  8 

WN  =  2(11  cjljfi  +  ||  ej| \  ]  =  2  [  cf  '  PXi  '  ci  +  '  ^r1  '  ei  ]  •  (3-9) 

t=l  i-l 

Because  WN  is  a  differentiable  function  of  Cj  and  e,-  for  all  non-zero  values  of  the 
residuals,  the  numerical  minimization  of  WN  subject  to  eq.  (3.2)  is  much  simpler  than 
the  minimization  of  W.  After  the  solution  is  obtained  one  may  check  the  membership 
values  fiXi  and  fiT  as  well  as  the  inequality  constraints  (3.4)  and  take  appropriate  action 
if  the  membership  values  are  not  positive  or  (3.4)  are  not  satisfied.  Usually  they  are 
satisfied  if  model  and  data  are  reasonable.  The  testing  of  the  conditions  is  most 
conveniently  done  in  terms  of  the  discords  D(Xi,XFi),  i.e.,  by  eq.  (3.5),  because  the 
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values  of  D  also  provide  a  measure  for  the  identification  of  outliers.  The  appropriate 
action  again  depends  on  the  particular  application  and  little  can  be  said  in  general.  It 
may  involve  a  removal  of  outliers,  reexamination  of  panderance  estimates,  change  of 
model,  etc. 

If  the  panderance  matrices  Pxi  and  PT  are  given,  then  eqs.  (3.2),  (3.4)  and  (3.9) 
define  a  constrained  minimization  problem  where  the  unknowns  are  T,  and  ef.  More 
common  are  situations  where  the  parameter  panderance  matrix  Pf,  of  the  free 
parameters  T  j  is  not  known.  Then  it  can  be  iteratively  determined  making  use  of  the 
panderance  propagation  formula  (2.41).  For  such  an  iteration,  one  assumes  an  initial 
estimate  Pjy0 ,  numerically  solves  the  minimization  problem  and  computes  the  ensuing 
panderance  matrix  PT,  x  using  the  linearized  panderance  propagation  formula  (2.41). 
Then  this  matrix  is  normalized  (see  below)  and  taken  as  initial  estimate,  etc. 

The  linearized  panderance  propagation  formula  produces  a  panderance  matrix 
with  a  norm  proportional  to  the  norms  of  the  input  panderances.  Therefore,  an 
iteration  without  a  normalization  of  the  iterated  panderance  matrix  between  iterations 
generally  does  not  converge  to  a  non-zero  finite  panderance  matrix  of  the  model 
parameters.  The  iteration  with  normalization  determines  on  the  other  hand  only  the 
relative  fuzziness  and  concordances  of  the  parameters  and  thereby  allows  one  to 
prescribe  the  relative  importance  of  model  and  data  fuzziness,  respectively.  Let  this 
relative  importance  be  expressed  by  a  scalar  parameter,  say  a.  Let  another  parameter, 
P,  express  the  relative  importance  of  the  fuzziness  of  the  free  and  prescribed  parameters, 
respectively.  We  express  the  objective  function  WN  as  a  function  of  the  two  parameters 
as  follows: 


s 


i=l 

(3.10) 

II 

O 

Z'** 

a 

)  P~l  ■  (C‘ 

V  le. 

(3.11) 

(1  -  a)2  ■  Pxi 

0 

0 

“  0  a2 

■PPti 

0 

(3.12) 

{  0  0  a2  ■  (l  —  0)2  ■  P Tp  j 

The  two  parameters  a  and  $  may  have  values  between  zero  and  one.  If  a=0  then 
we  have  an  adjustment  problem  with  fuzzy  data  and  crisp  model.  If  a  =  l  then  the  data 
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are  crisp  and  the  model  is  fuzzy.  If  a>0  and  $=0  then  the  data  are  fuzzy,  the  free 
model  parameters  are  crisp  and  the  prescribed  model  parameters  are  fuzzy.  If  or>0  and 
0—1  then  the  data  are  fuzzy  and  and  the  model  has  only  free  and  fuzzy  parameters.  By 
choosing  appropriate  values  for  a  and  0  one  can  obtain  a  number  of  different 
combinations  of  conditions.  We  give  three  examples  of  such  combinations.  First,  if  all 
the  panderance  matrices  P Xi<  P Tf  and  PTp  are  prescribed,  then  we  use  (3.12)  with 
Qr=2/3  and  $=1/2.  Second,  if  PXi  and  PTp  are  given,  and  PTj  is  iteratively  determined, 
then  we  specify  a  and  0  by  two  conditions.  First,  we  require  that  the  relative  weights  of 
the  contributions  Ell  C.\\h  and  EIMS-,  to  WN  are  the  same.  Second,  we  require 
that  the  contribution  of  Ell  efi  II  \f  (after  the  normalization  of  the  iterated  parameter 
panderance  matrix  PTj )  to  have  the  same  weight  as  the  sum  of  the  contributions  of 
Elklljw  and  ElUpilliy  The  first  requirement  is  satisfied  if 


a  = 


1 

2-0  ‘ 


The  second  requirement  we  express  by  the  equation 


(3.13) 


a  a  a 

E  IkllJ,  +  E  II e, ill f,  -  E  IM?,  -0  .  (3.14) 

i=l  i=l  i'=l 

Eqs.  (3.13)  and  (3.14)  are  a  system  of  equations  for  a  and  0,  because  the  residuals 
depend  on  a  and  0.  We  can  obtain  the  numerical  solution  of  this  system,  for  instance, 
using  a  regula  falsi  algorithm  for  0. 


As  a  third  example  for  the  choice  of  a  and  0  we  assume  that  the  PXi  are  known 
and  PTj  is  iteratively  determined,  as  before,  and  that  PTp  is  known  only  up  to  scaling 
factor.  Now,  PTj  as  well  as  PTp  contain  arbitrary  scaling  factors.  In  this  case,  we 
postulate  the  following  two  conditions  for  a  and  0.  First,  we  require  that  the 
contribution  of  Ell e/, Hr/  to  kv  receives  the  same  weight  as  the  contribution  of 
EIMV  Second,  we  require  that  the  sum  of  all  contributions  of  the  parameter 
residuals  has  the  same  weight  as  the  contribution  Ell  c,l \h  of  the  residuals  of  the 
observations.  The  two  conditions  can  be  expressed  in  form  of  the  following  two 
equations. 


and 


C(M)-S  II^IISy-E  lie, ill!-,  =0 

1=1  1=1 


(3.15) 
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(3-16) 


«(»,«  — E  wi+S  1.,^+E MK-o  . 


i=l 


i=l 


«=1 


This  system  of  equations  for  or  and  0  we  solve  numerically  by  a  regula  falsi  algorithm  in 
the  a,/#-plane.  (Note  that  G(a, 0)<0,  G(a,  1)>0,  H(O,0)<O  and  H(1,0)>O,  which  assures 
the  existence  of  a  solution  within  the  unit  square). 


3.3.  Numerical  Treatment. 


The  central  part  of  the  numerical  treatment  of  the  model  fitting  problem  is  the 
solution  of  the  following  constrained  minimization  problem: 


Minimize 


subject  to 


W^a.0)  =  E  («f,  ef)  />-■  *<  , 

1=  1  l.C»J 

F i(Xi  +  c,- ,  T  +  e,)  =  0  ,  i'  =  l . s  , 


(3.17a) 

(3.17b) 


where  the  Pfi  are  defined  by  eq.  (3.12).  This  problem  formally  is  equal  to  a  usual  least 
squares  model  fitting  problem  with  vector  data  and  implicit  model  equations  [2j. 
Software  for  the  solution  of  such  problems  is  available  for  instance  from  [l]  and  can  be 
used  to  solve  (3.17).  The  solution  of  the  example  presented  in  Section  4  was  computed 
using  the  utility  routine  COLSAC  [1],  The  result  of  a  COLSAC  calculation  consists  of 
the  optimal  values  of  cit  «,•  and  T,  and  a  posterior  estimate  of  the  panderance  matrix 
PTf  computed  with  the  panderance  propagation  formula  (2.41)  from  the  input 
panderance  matrices.  If  no  iterations  are  involved,  i.e.,  if  PT  is  prescribed,  then  a  single 
call  of  COLSAC  concludes  the  calculations. 


The  inequality  constraints  (3.4)  usually  are  inactive,  except  for  outliers.  If  it  is 
determined  that  data  with  large  discords  cannot  be  discarded  or  treated  in  the  same 
manner  as  other  data,  (i.e.,  the  inequality  constraints  ignored),  then  one  has  a  more 
complicated  programming  problem  at  hand.  There  is  no  guarantee  that  the  problem 
has  a  solution  and,  if  a  solution  exists,  it  generally  is  not  very  useful,  because  it 
particularly  accommodates  outliers.  The  treatment  of  problems  with  outliers  and  too 
stringent  conditions  (3.4)  very  much  depends  on  the  application  and  typically  may 
involve  a  change  of  the  model  functions  F We  shall  not  further  pursue  this  type  of 
problem,  and  assume  that  the  inequality  constraints  are  inactive  at  the  solution. 

For  certain  ranges  of  a  and  0  the  convergence  of  the  procedure  for  the 
determination  of  PT  can  be  quite  slow.  In  such  cases,  one  may  first  determine  the 
matrix  PT  by  iteration  for  extreme  values  of  a  and  0  and  use  an  interpolated  matrix  PT 


-24  - 


for  subsequent  calculations  with  other  a  and  /?.  Then,  iteration  to  find  the  final  PT 
might  be  resumed  in  the  vicinity  of  the  final  point  in  the  a, /3-plane. 


4.  MODEL  OF  PLATE  DEFORMATION. 


4.1  Data  and  Model. 

The  vulnerability  of  armored  vehicles  depends  among  other  things  on  the 
deformation  of  armor  plates  by  the  blast  of  a  nearby  explosion.  References  [10]  and  [11] 
describe  a  series  of  experiments  in  which  such  deformations  were  measured.  We  shall 
use  data  from  these  experiments  in  an  example  of  fuzzy  model  fitting. 

The  relevant  data  consist  of  the  charge  energy,  the  distance  of  the  charge  from  the 
armor  plate,  the  thickness  of  the  plate  and  the  maximum  deflection  of  the  plate.  (The 
size  of  the  unsupported  span  of  the  plate  was  the  same  in  all  experiments  and,  therefore, 
its  influence  on  the  deformation  cannot  be  obtained  from  the  data.  However,  in  most  of 
the  observed  cases,  the  size  of  the  deformation  was  such  that  the  influence  of  the  plate’s 
support  reasonably  could  be  assumed  negligible).  Some  of  these  data  were  measured 
quite  precisely  (e.g.,  the  plate  thickness),  while  other  data  were  inaccurate  and  obvious 
candidates  for  fuzzy  numbers.  Such  a  datum  is,  e.g.,  the  strength  of  the  plate’s  material 
which  we  chose  to  characterize  by  its  ultimate  yield  stress.  The  yield  stress  is  ill  defined 
and  difficult  to  measure,  and  it  was  found  to  vary  considerably  between  plate  specimens 
[llj.  In  summary,  it  was  assumed  that  all  data  were  fuzzy  numbers  with  appropriate 
spreads. 

The  regression  model  of  the  plate  deformation  is  not  derived  from  a  mechanical 
theory  of  the  event,  because  such  a  general  theory  does  not  exist.  Instead,  the  model 
was  determined  by  exploratory  data  analysis  and  dimensional  arguments.  This 
uncertainty  about  the  form  of  the  model  was  represented  in  the  model  specification  by 
using  fuzzy  equations  as  model  constraints.  Thus,  in  this  example,  the  model  as  well  as 
the  data  are  fuzzy. 


The  following  two  dimensionless  quantities  were  used  as  variables: 


A 


B  = 


E 

d2  Tov 


(4.1) 


where  z  (m)  is  the  maximum  deflection,  d  (m)  is  the  distance  of  the  charge  from  the 
plate,  E  (J)  is  the  energy  released  by  the  explosion,  T  (m)  is  the  thickness  of  the  plate 
and  <rt  (Pa)  is  the  ultimate  yield  stress  of  the  plate’s  material.  The  maximum  deflection 
was  obtained  from  a  detailed  analysis  of  various  observed  deflections,  as  described  in 
reference  [5],  d,  E  and  T  were  obtained  by  simple  measurements,  and  for  <ry  a  nominal 
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handbook  value  was  used. 

In  reference  [5]  the  quantity  B  was  called  a  ballistic  limit  indicator,  because  the 
size  of  B  is  a  predictor  for  the  onset  of  plate  failure:  the  plate  fails  (develops  a  hole)  if  B 
exceeds  the  fuzzy  bound  £0={129.,40.2}.  Since  B  is  also  useful  as  an  indicator  for  the 
size  of  plate  deformation  by  blast,  a  better  name  for  B  is  the  blast  damage  indicator. 

Table  1  lists  the  deflection  data  with  their  estimated  spreads  and  concordances. 
The  spreads  of  d,  E,  T  and  av  were  assumed  to  be  constant  and  estimated  by 
reasonable  assumptions  about  measurement  accuracy.  The  spread  of  the  deflection  z 
was  obtained  by  the  data  analysis  in  reference  [5].  The  spreads  of  A  and  B  were 
computed  from  these  estimates  using  eq.  (2.41).  The  concordances  between  the  observed 
^  and  the  corresponding  B  are  not  zero,  because  both  depend  on  the  distance  d. 
Figures  2  and  3  display  the  data  points  without  and  with  the  corresponding  supports  of 
their  membership  functions,  respectively.  The  different  symbols  of  data  points  in  the 
figures  distinguish  different  subsets  of  the  experiments  but  the  distinction  is  not 
important  for  the  present  example. 

The  model  function  was  chosen  to  be 

F{X,B  ]  a,b,c)  =  a  +  b  ■  B*  -  A  .  (4.2) 

The  parameters  a  and  b  were  assumed  to  be  free  and  their  spreads  were  not  prescribed. 
The  parameter  <T  was  assumed  to  have  the  apex  1.50  and  either  a  spread  of  0.2  or  an 
undetermined  spread.  Thus,  in  this  example,  the  panderance  matrix  PTJ  of  the  free 
parameters  is  a  (2X2)-matrix  and  the  panderance  matrix  PTp  of  the  prescribed 
parameter  c  is  a  scalar.  The  data  fitting  was  carried  out  as  described  in  Section  3  with 
the  aid  of  the  utility  routine  COLSAC  [lj. 


4.2  Data  and  Model  Spreads. 

The  data  panderance  matrices  PXi  are  given  in  terms  of  the  data  in  Table  1  by 

sAi  5/ti  '  sBi  '  cABi 

sAi  '  sBi  '  cABi 


PXi  = 


(4.3) 


where  and  sBi  are  the  spreads  of  A,  and  Bit  respectively  and  cABi  is  the  concordance 
between  A,  and  B{.  The  parameter  part  of  the  regression  panderance  matrix  (3.12)  is 


Pr  = 


a2  ■  0®  •  PT; 

0 


0 

a2  (l-  pf  ■  0.22 


(4.4) 
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where  the  (2X2)-matrix  PTj  was  determined  by  iteration.  We  calculated  two  cases.  In 
one  case,  PT  was  assumed  as  given  by  eq.  (4.4).  In  the  other  case,  we  assumed  that  the 
spread  of  the  exponent  c  (  =0.2  in  eq.  (4.4) )  is  not  known,  that  is,  PTp  was  assumed  to 
contain  an  unknown  scaling  factor. 

In  the  former  case,  the  normalization  factor  of  PTj  was  chosen  such  that  the 
maximum  discord  between  any  of  the  apexes  AT,-  of  the  observations  X{,  i  =  l,...,s,  and 
the  fuzzy  surface  F(X,  T  j,  Tp)=0  was  equal  to  five.  In  the  latter  case,  we  first 
normalized  PTj  and  PTp  separately  such  that  the  maximum  discords  between  the 
observation  apexes  and  F(X,T  /,Tp)= 0  and  F(X,T f,T  p)  = 0,  respectively,  were  equal  to 
five.  Then  the  complete  matrix  PT  (consisting  of  the  two  separately  normalized 
submatrices)  was  normalized  such  that  the  maximum  discord  between  the  observation 
apexes  and  the  surface  F(X,T  j,  Tp)=0  was  equal  to  five.  We  remind  that  the  choice  of 
a  normalization  procedure  (in  the  present  case,  also  the  maximum  discord  values  chosen 
for  the  normalization)  is  arbitrary  in  the  sense  that  it  affects  only  the  numerical  values 
of  a  and  /?,  which  regulate  the  relative  weights  of  data  and  parameter  residuals  (see 
Section  3.2). 

The  posterior  fuzziness  of  the  fitted  model  can  be  estimated  using  the  panderance 
propagation  formula  (2.41)  if  the  panderances  of  data  and  model  parameters  are  all 
prescribed.  If  a  part  of  the  model  parameter  panderance  matrix  PT  is  iteratively 
determined  then  the  result  contains  arbitrary  parameters  and  a  direct  application  of 
(2.41)  cannot  be  justified.  Therefore,  one  should  in  such  cases  determine  the  spread  of 
the  model  by  other  means,  also  taking  into  account  the  actually  observed  deviations 
between  data  and  fitted  model. 

In  the  present  example,  the  fitted  function  was  used  in  the  form 

A{B)=A{B,T)  =  a+  b-  B*  ,  (4.5) 

and  its  spread  was  adapted  to  observations  by  using  two  quadratic  modulator  functions 
M(B),  one  each  for  the  positive  and  negative  spreads  of  A.  First,  PT  was  normalized 
using  a  factor  V  such  that  the  largest  discord  between  the  fitted  curve  A(B,T)  and  the 
Xi  was  five.  This  gave  a  reasonable  spread  increase  pattern  for  the  model  outside  the 
observed  range  of  B  but  a  too  optimistic  spread  within  the  observed  range.  Therefore, 
modulator  functions  M(B)  were  introduced  to  define  modulated  panderance  matrices  by 

=  -  1*  PT  ■  (4.6) 

The  modulator  functions  were  determined  such  that  the  largest  discord  between  the 

«■*»»  ~ 

fitted  A(B,T)  and  the  observations  AT,-  was  less  than  two.  The  form  of  the  modulator 
functions  was 
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M{B)  =  max  {  1 ,  m1  +  m2  •  (B  —  J50)2  } 


(4.7) 


with  proper  constants  Ttiv  m2  and  B0,  and  a  different  function  was  determined  for 
positive  and  negative  differences,  respectively,  between  the  fitted  function  and 
observations.  The  maximal  discord  value  of  two  was  established  by  numerical 
experiments  and  chosen  because  it  produced  a  spread  pattern  in  reasonable  agreement 
with  the  scatter  of  observations.  The  final  formula  for  the  spread  of  A(B)=A(B,  T)  is 


Sa  = 


dA  r 

dA 

T 

^  *T  mod 

dT 

1/2 


(4.8) 


4.3.  Numerical  Results. 

Table  2  contains  the  numerical  results  in  the  two  presented  examples.  The  table 
lists  the  optimal  values  of  the  model  parameters,  their  spreads  and  the  concordance  cab 
between  a  and  6.  The  panderance  matrix  of  the  parameter  vector  is  given  in  terms  of 
the  spreads  and  cab  by  eq.  (2.6).  The  quoted  spreads  correspond  to  a  matrix  that  is 
normalized  such  that  in  the  parameter  space  the  largest  discord  between  the  locally 
corrected  parameters  and  the  apex  of  the  parameter  vector  equals  one  and,  therefore, 
the  listed  spread  sc  of  c  in  Case  1  is  not  equal  to  the  prescribed  nominal  value  of  0.2  .  A 
comparison  of  the  two  cases  illustrates  the  effect  of  assigning  a  fixed  value  to  the  spread 
of  a  model  parameter.  Thus,  in  the  present  example,  the  spread  of  c  is  larger  in  Case  1 
than  in  Case  2,  as  are  the  uncertainties  of  the  other  two  model  parameters  a  and  b.  The 
increase  of  the  latter  two  is  a  consequence  of  the  condition  (3.16).  Plots  of  the  final 
model  curve  X(B)  are  shown  in  Figures  4  and  5.  The  figures  display  for  each  data  point 
its  corresponding  support  ellipse  (modified  by  the  factor  1— a,  see  eq.  (3.12)),  the 
corrected  data  point  X^+c,,  and  a  short  segment  of  the  fitted  function  that  passes 
through  the  corrected  point  and  has  the  highest  membership  value  nFT.  (It  is  a  segment 
of  the  curve  F(X,  T  +  e,)=0  ).  The  spreads  of  the  model  curves  shown  in  Figures  4  and  5 
are  computed  using  the  factor  if  and  the  modulator  functions  M(B)  from  Table  2  to 
modulate  the  parameter  panderance  matrix  given  at  the  top  of  the  table  (see 
Section  4.2).  The  differences  between  the  shapes  of  the  spreads  of  the  fitted  curves  in 
both  figures  illustrate  how  a  prescribed  model  fuzziness  may  affect  the  result  of  the 
fitting. 

An  overview  of  the  local  parameter  adjustments  is  provided  by  Figures  6  and  7. 
They  show  in  the  parameter  space  three  views  of  the  locally  adjusted  parameter  sets, 
i.e.,  the  end  points  of  the  vectors  r  +  e,=(o  +  e4t-  ,b+ebi  ,c+eM)r.  Also  shown  are  the 
contours  of  the  support  ellipsoid  of  the  fuzzy  parameter  vector  T.  The  ellipsoid 
corresponds  to  the  spreads  sa,  sb  and  ae  listed  in  Table  2. 
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5.  SUMMARY. 

Mathematical  models  of  terminal  ballistics  problems  often  must  be  based  on 
incomplete  theories  and  data  that  contain  inaccuracies  which  are  not  probabilistic.  An 
example  of  such  a  problem  is  the  deformation  of  an  armor  plate  by  a  nearby  explosion. 
To  treat  this  problem,  we  have  used  concepts  of  fuzzy  set  theory  and  developed 
algorithms  for  a  special  type  of  fuzzy  regression. 

The  algorithms  are  based  on  concepts  of  analytic  geometry  and  restricted  to  data 
and  model  parameter  vectors  which  are  characterized  by  a  particular  type  of  conical 
membership  function.  The  regression  principle  is  to  maximize  in  a  least  squares  sense 
membership  values  of  fuzzy  observations  and  of  fuzzy  fitting  functions.  The  restriction 
to  conical  membership  functions  and  the  employment  of  a  least  squares  principle  have 
the  advantage  that  the  calculations  can  be  done  with  available  software  for  least  squares 
regression  problems  with  implicit  constraint  functions.  The  model  functions  can  be 
implicit  and  nonlinear  with  respect  to  data  as  well  as  with  respect  to  model  parameters, 
and  the  algorithms  accept  fuzzy  data  as  well  as  fuzzy  parameters,  whereby  the  latter 
may  be  either  free  or  prescribed.  The  method  is  more  flexible  than  previously  published 
methods,  which  are  restricted  either  to  linear  models,  [6]  through  [9],  [12]  and  [13],  or  to 
crisp  models  [3],  or  to  crisp  data  [4[.  (The  membership  functions  are  restricted  to 
particular  types  in  all  of  the  quoted  publications  about  fuzzy  regression  methods).  We 
feel  that  the  possibility  to  choose  freely  the  model  functions  is  a  considerable  asset  of  the 
present  method. 

The  method  is  applied  to  a  set  of  data  from  experiments  where  armor  plates  were 
deformed  by  bare  charge  explosions.  The  result  is  a  fuzzy  non-linear  function  that 
predicts  the  maximum  of  a  relative  deflection  of  the  plate  in  terms  of  a  dimensionless 
blast  damage  indicator.  The  model  takes  into  account  the  inaccuracies  of  data  and  the 
vagueness  of  the  theoretical  basis  by  assigning  proper  fuzziness  measures  to  the 
prediction. 
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e  2.  Observed  plate  deflection  versus  load 
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Figure  6.  Local  parameter  adjustments  in  Case  1 
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Figure  7.  Local  parameter  adjustments  in  Case  2 


Table  1.  Input  Data, 


Nr. 

Label 

B 

B-Spread 

A 

A- Spread 

Concordance 

1 

RD-  2/T#  5 

105.26 

21.17 

0.76056 

001369 

0  053516 

MM 

RD-  5/T#  5 

102.93 

20.70 

0.76508 

0.01158 

0.063669 

RD-  6/T#  5 

11821 

23.81 

0.98405 

0.01693 

0.063657 

wm 

RD-  8/T#  5 

51.00 

10  21 

0.43281 

0.00458 

0.018055 

5 

RD-  9/T#  5 

64.24 

12.87 

0.55264 

0.00386 

0.034592 

6 

RD-10/T#  5 

51.34 

10.28 

0.37383 

0.00344 

0.034592 

7 

RD-ll/T#  5 

64.99 

43.02 

0  41531 

0.00366 

0.029891 

8 

RD-12/T#  5 

65.18 

13.06 

0.43272 

0  00582 

0001796 

9 

RD-13/T#  5 

85  16 

17.07 

0.51172 

0.00408 

0.039626 

10 

RD-14/T#  5 

115.30 

23.12 

086140 

0  00971 

0.038086 

11 

RD-16/T#  5 

123.73 

24  82 

1.02472 

0.02159 

0.022186 

12 

RD-18/T#  5 

97  93 

19  63 

0.77744 

0.00821 

0.034652 

13 

RD-19/T#  5 

107  91 

21  64 

0  87198 

0.00772 

0044735 

14 

RD-21/T#  5 

81.80 

16  39 

0  67031 

0.00425 

0  030151 

15 

RD-22/T#  5 

87.36 

17.50 

0.77164 

0.00593 

0  026309 

16 

RD-23/T#  5 

89  61 

17.95 

0.79166 

000471 

0  036016 

17 

RD-24/T#  5 

101  66 

20.37 

0  92772 

0  00577 

0  038845 

18 

RD-25/T#  5 

112  61 

22.56 

0.91424 

0  01386 

0016990 

19 

RD-27/T#  5 

119.06 

23  85 

092703 

0  00606 

0  039406 

20 

RD-  1/T#10 

69  83 

14  07 

0.40709 

0.00767 

0  066639 

21 

RD-  2/T#  10 

135  94 

27.54 

0.96792 

0.02033 

0  116586 

22 

RD-  3/T#10 

95  37 

19  25 

0.62306 

001206 

0  088135 

23 

RD-  4/T#10 

116  19 

23  51 

0  77957 

0.01784 

0  088517 

24 

RD-  7/T#10 

96  04 

19  39 

0  69616 

0.01057 

0  112396 

25 

RD-  8/T#  10 

114  02 

23  05 

0  78626 

0.01431 

0.111430 

26 

RD-  9/T#  10 

114  04 

23  06 

078291 

001279 

0  124088 

27 

RD-10/T#10 

131  31 

2670 

1  28464 

0  03112 

0  101073 

28 

RD-11/T#10 

129  67 

26.07 

0  95972 

0  01227 

0  075380 

29 

RD-14/T#10 

122  87 

24.70 

1.02897 

001251 

0  079290 

30 

RD-16/T#10 

130  02 

26  09 

1.19877 

0.01199 

0061773 

31 

RD-18/T#10 

130  85 

26.26 

1.13072 

001138 

0  061376 

32 

RD-19/T#10 

152.65 

30  64 

1.16320 

0  01258 

0  063302 

33 

RD-22/T#10 

123  06 

24  68 

1.05742 

001652 

0.032715 

34 

RD-  1/T#20 

86  41 

17  55 

0.20707 

0  00873 

0.071527 

35 

RD-  2/T#20 

11083 

22  60 

0.63734 

0.01885 

0  129906 

36 

RD-  3/T#20 

66  25 

13.32 

0.28792 

0  00438 

0  063388 

37 

RD-  4/T#20 

86.75 

17  46 

0.63775 

0.00852 

0094093 

38 

RD-  6/T#20 

100.05 

20.16 

0  54149 

0.00845 

0.093388 

39 

RD-  8/T#20 

115  95 

23  30 

0.65883 

001016 

0.062502 

40 

RD-12/T#20 

90  62 

18  19 

0  45640 

0.00668 

0  052095 

41 

RD-13/T#20 

102  40 

20  57 

0  59983 

0  00682 

0  075114 

42 

RD-15/T#20 

91  35 

18  34 

0.48248 

0.00518 

0  071063 

43 

RD-16/T#20 

90  17 

18  15 

0  40348 

0  00625 

0.081118 

44 

RD-17/T#20 

122  35 

24  69 

0  60291 

0  00953 

0  107990 

Table  2.  Model  Parameters. 

Case  1.  Prescribed  Apex  and  Spread  of  <?. 

a  =  {  —  5.9210-2,  31.6010-2},  sfl -^  =  5.83  10"2 

$={  7.52  10~4,  3.68  10-4  },  sb =  0.68- 10-4 

c  =  {  1.500  ,  0.134  }  ,  V  ^  =  0.025 

cab  =  -  0.952  52,  *  =  0.184  39 

M(B)+  =  2.588  -  1.19-10~4  (B  -  48.4)2 

M{B)_  =  -  2.617  +  3.75  TO ~4  (B  -  90.8)2 

a  =  0.75,  0  =  0.67 

Case  2.  Prescribed  Apex  and  Undetermined  Spread  of  "c. 

a  =  {  —  3.01  T0~2,  20.88  TO-2  }  ,  ^  =  12.28  10~2 

£={  7.50  10~4,  2.28  10“4  },  sb-  'if  —  1.5210-4 

c  =  {  1.500  ,  0.040  }  ,  v  $  =  0.024 

cab  =  -0.955  21,  ^  —  0.587  88 

M(B)+  =  1.952  -  4.18  10_4(B  -  80.4)2 

M{B)_  =  -  2.641  +  3.80  10 ~4  {B  -  93.4)2 

a  =  0.60,  0  =  0.43 

Note:  The  fuzzy  numbers  in  this  table  are  given  in  the  form  {a,sa},  that  is,  the 
second  number  in  the  braces  is  the  spread  (not  the  panderance). 
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No  of  No  of 

Copies  Organization  Copies  Organization 


1  Office  of  the  Secretary  of  Defense 
OUSD(A) 

Director,  Live  Fire  Testing 
ATTN:  James  F.  O’Bryon 
Washington,  DC  20301-3110 

2  Administrator 

Defense  Technical  Info  Center 
ATTN:  DTIC-DDA 
Cameron  Station 
Alexandria,  VA  22304-6145 

1  HQDA  (SARD-TR) 

WASH  DC  20310-0001 

1  Commander 

US  Army  Materiel  Command 
ATTN:  AMCDRA-ST 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333-0001 

1  Commander 

US  Army  Laboratory  Command 
ATTN:  AMSLC-DL 
Adelphi,  MD  20783-1145 

2  Commander 

US  Army,  ARDEC 

ATTN:  SMCAR-IMI-I 

Pica  tinny  Arsenal,  NJ  07806-5000 

2  Commander 

US  Army,  ARDEC 

ATTN:  SMCAR-TDC 

Pica  tinny  Arsenal,  NJ  07806-5000 

1  Director 

Benet  Weapons  Laboratory 
US  Army,  ARDEC 
ATTN:  SMCAR-CCB-TL 
Watervliet,  NY  12189-4050 

1  Commander 

US  Army  Armament,  Munitions 
and  Chemical  Command 
ATTN:  SMCAR-ESP-L 
Rock  Island,  IL  61299-5000 

1  Commander 

US  Army  Aviation  Systems  Command 
ATTN:  AMSAV-DACL 
4300  Goodfellow  Blvd. 

St.  Louis,  MO  63120-1798 


1  Director 

US  Army  Aviation  Research 
and  Technology  Activity 
ATTN:  SAVRT-R  (Library) 

M/S  219-3 

Ames  Research  Center 
Moffett  Field,  CA  94035-1000 

1  Commander 

US  Army  Missile  Command 
ATTN:  AMSMI-RD-CS-R  (DOC) 
Redstone  Arsenal,  AL  35898-5010 

1  Commander 

US  Army  Tank-Automotive  Command 
ATTN:  AMSTA-TSL  (Technical  Library) 
Warren,  MI  48397-5000 


1  Director 

US  Army  TRADOC  Analysis  Command 
ATTN:  ATAA-SL 

White  Sands  Missile  Range,  NM  88002-5502 


(ciau.  only)  j  Commandant 

US  Army  Infantry  School 
ATTN:  ATSH-CD  (Security  Mgr.) 
Fort  Benning,  GA  31905-5660 

(Uaduo.  **'»)  1  Commandant 

US  Army  Infantry  School 
ATTN:  ATSH-CD-CSO-OR 
Fort  Benning,  GA  31905-5660 


1  Air  Force  Armament  Laboratory 
ATTN:  AFATL/DLODL 
Eglin  AFB,  FL  32542-5000 


Aberdeen  Proving  Ground 

2  Dir,  USAMSAA 

ATTN:  AMXSY-D 

AMXSY-MP,  H.  Cohen 

1  Cdr,  USATECOM 

ATTN:  AMSTE-TD 

3  Cdr,  CRDEC,  AMCCOM 

ATTN:  SMCCR-RSP-A 
SMCCR-MU 
SMCCR-MSI 

1  Dir,  VLAMO 

ATTN:  AMSLC-VL-D 
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Commander 

US  Army  Laboratory  Command 
ATTN:  AMSLC-CG 

AMSLC-TR,  R.  Vitali 
2800  Powder  Mill  Road 
Adelphi,  MD  20783-1145 

Commander 

US  Army  Missile  Command 
ATTN:  AMSTA-CG 
Redstone  Arsenal,  AL  35898-5000 

Commander 

US  Army  Concepts  Analysis  Agency 
ATTN:  D.  Hardison 
8120  Woodmont  Avenue 
Bethesda,  MD  20014-2797 

C.I.A. 

10  lR/DB/Standard 
GE47  HQ 

Washington,  DC  20505 

Commander 
US  Army  War  College 
ATTN:  Library-FF229 
Carlisle  Barracks,  PA  17013 

US  Army  Ballistic  Missile 
Defense  Systems  Command 
Advanced  Technology  Center 
P.O.  Box  1500 
Huntsville,  AL  35807-3801 

Commander 

US  Army  Materiel  Command 
ATTN:  AMCPM-GCM-WF 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333-5001 

Commander 

US  Army  Materiel  Command 
ATTN:  AMCDE-DW 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333-5001 


1  Commander 

US  Army  Watervliet  Arsenal 
ATTN:  SARWV-RD,  R.  Thierry 
Watervliet,  NY  12189-5001 

3  Commander 

US  Army  Natick  R&D  Center 
ATTN:  STRNC-OI 
Natick,  MA  01760 

1  Commander 

US  Army  Aviation  Systems  Command 
ATTN:  AMSAV-ES 
4300  Goodfellow  Blvd. 

St.  Louis,  MO  63120-1798 

1  Commander 

US  Army  Survivability  Management  Office 
ATTN:  SLCSM-D 
2800  Powder  Mill  Road 
Adelphi,  MD  20783-1145 

1  Commander,  USACECOM 

R&D  Technical  Library 
ATTN:  ASQNC-ELC-I-T,  Myer  Center 
Fort  Monmouth,  NJ  07703-5000 

1  Commander 

US  Army  Harry  Diamond  Laboratories 
ATIN:  SLCHD-TA-L 
2800  Powder  Mill  Road 
Adelphi,  MD  20783-1145 

1  Commandant 

US  Army  Aviation  School 
ATTN:  Aviation  Agency 
Fort  Rucker,  AL  36360 

1  Project  Manager 

US  Army  Tank-Automotive  Command 
Improved  TOW  Vehicle 
ATTN:  AMCPM-ITV 
Warren,  MI  48397-5000 

2  Project  Manager 

Ml  Abrams  Tank  System 

ATTN:  AMCPM-ABMS-SA,  T.  Dean 

Warren,  MI  48092-2498 


No.  of 

Copies  Organization 

1  Project  Manager 

Fighting  Vehicle  Systems 
ATTN:  AMCPM-BFVS 
Warren,  MI  48092-2498 

1  President 

US  Army  Armor  &  Engineer  Board 
ATTN:  ATZK-AD-S 
Fort  Knox,  KY  40121-5200 

1  Commander 

US  Army  Vulnerability  Assessment 
Laboratory 
ATTN:  SLCVA-CF 

White  Sands  Missile  Range,  NM  88002-5513 

1  Project  Manager 

M-60  Tank  Development 
ATTN:  AMCPM-ABMS 
Warren,  MI  48092-2498 

1  Commander 

US  Army  Training  &  Doctrine  Command 
ATTN:  ATCD-MA,  MAJ  Williams 
Fort  Monroe,  VA  23651 

2  Commander 

US  Army  Materials  Technology  Laboratory 
ATTN:  SLCMT-ATL 
Watertown,  MA  02172-0001 

1  Commander 

US  Army  Research  Office 
ATTN:  Technical  Library 
P.O.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 

1  Commander 

US  Army  Belvoir  R&D  Center 

ATTN:  STRBE-WC 

Fort  Belvoir,  VA  22060-5606 

1  Director 

US  Army  TRAC  -  Ft.  Lee 
Defense  Logistics  Studies 
Fort  Lee,  VA  23801-6140 


No.  of 

Copies  Organization 
1  Commandant 

US  Army  Command  and  General  Staff 
College 

Fort  Leavenworth,  KS  66027 

1  Commandant 

US  Army  Special  Warfare  School 
ATTN:  Rev  &  Tng  Lit  Div 
Fort  Bragg,  NC  28307 

1  Commander 

US  Army  Foreign  Science  &  Technology 
Center 

ATTN:  AMXST-MC-3 
220  Seventh  Street,  NE 
Charlottesville,  VA  22901-5396 

2  Commandant 

US  Army  Field  Artillery  Center  &  School 
ATTN:  ATSF-CO-MW,  B.  Willis 
Fort  Sill,  OK  73503-5600 

1  Office  of  Naval  Research 
ATTN:  Code  473,  R.S.  Miller 
800  N.  Quincy  Street 
Arlington,  VA  22217-9999 

3  Commandant 

US  Army  Armor  School 

ATTN:  ATZK-CD-MS,  M.  Falkovitch 

Armor  Agency 

Fort  Knox,  KY  40121-5215 

2  Commander 

Naval  Sea  Systems  Command 
ATTN:  SEA62R 
SEA  64 

Washington,  DC  20362-5101 

1  Commander 

Naval  Air  Systems  Command 
ATTN:  AIR-954-Tcchnical  Library 
Washington,  DC  20360 

1  Naval  Research  Laboratory 
Technical  Library 
Washington,  DC  20375 
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No.  of 

Copies  Organization 

1  Commander 

Naval  Surface  Warfare  Center 
ATTN:  Code  DX-21  Technical  Library 
Dahlgren,  VA  22448-5000 

1  Commander 

Naval  Weapons  Center 

ATTN:  Information  Science  Division 

China  Lake,  CA  93555-6001 

1  Commander 

Naval  Ordnance  Station 
ATTN:  Technical  Library 
Indian  Head,  MD  20640-5000 

1  AF  Astronautics  Laboratory 

AFAL/TSTL,  Technical  Library 
Edwards  AFB,  CA  93523-5000 

1  AFATL/DLYV 

Eglin  AFB,  FL  32542-5000 

1  AFATUDLXP 

Eglin  AFB.  FL  32542-5000 

1  afativdue 

Eglin  AFB,  FL  32542-5000 

1  NASA/Lyndon  B.  Johnson  Space  Center 
ATTN:  NHS22  Library  Section 
Houston,  TX  77054 

1  FTD/NBS 

Wright-Patterson  AFB,  OH  45433 

1  Director 

Lawrence  Livermore  Laboratory 
ATTN:  Technical  Information 
Department  L-3 
P.O.  Box  808 
Livermore,  CA  94550 

2  Director 

Los  Alamos  Scientific  Laboratory 
ATTN:  Document  Control  for  Reports 
Library 

P.O.  Box  1663 

Los  Alamos,  NM  87544 


No.  of 

Copies  Organization 

1  Director 

Sandia  Laboratories 
ATTN:  Document  Control  for  3141 
Sandia  Report  Collection 
Albuquerque,  NM  87115 

1  Director 

Sandia  Laboratories 
Livermore  Laboratory 
ATTN:  Document  Control  for  Technical 
Library 
P.O.  Box  969 
Livermore,  CA  94550 

1  Director 

National  Aeronautics  and  Space 
Administration 

Scientific  &  Technical  Information 
Facility 
P.O.  Box  8757 

Baltimore/Washington  International  Airport, 
MD  21240 

1  Aerospace  Corporation 

ATTN:  Technical  Information  Services 

P.O.  Box  92957 

Los  Angeles,  CA  90009 

1  The  Boeing  Company 

ATTN:  Aerospace  Library 
P.O.  Box  3707 
Seattle,  WA  98124 

1  Director 

Institute  for  Defense  Analyses 
ATTN:  Library 
1801  Beauregard  Street 
Alexandria,  VA  22311 
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Battelle  Memorial  Institude 
ATTN:  Technical  Library 
505  King  Avenue 
Columbus,  OH  43201-2693 

Johns  Hopkins  University 
Applied  Physics  Laboratory 
ATTN:  Jonathan  Fluss 
John  Hopkins  Road 
Lauel,  MD  20707-0690 

Pennsylvania  State  University 
Department  of  Mechanical  Engineering 
ATTN:  K.  Kuo 
University  Park,  PA  16802-7501 

SRI  International 
Propulsion  Sciences  Division 
ATTN:  Technical  Library 
333  Ravenswood  Avenue 
Menlo  Park,  CA  94025-3493 

Rensselaer  Polytechnic  Institute 
Department  of  Mathematics 
Troy,  NJ  12181 

Eichelberger  Consulting  Company 
ATTN:  Dr.  R.  Eichelberger,  President 
409  West  Catherine  Street 
Bel  Air,  MD  21014 

FMC  Corporation 
Aberdeen  Regional  Office 
ATTN:  Leland  Watermeier 
Bel  Air,  MD  21014 


Aberdeen  Proving  Ground 

12  Dir,  USAMSAA 

ATTN:  AMXSY-A,  W.  Clifford 
J.  Meredith 
AMXSY-C,  A.  Reid 
AMXSY-CS,  P.  Beavers 

C.  Cairns 

D.  Frederick 
AMXSY-G,  J.  Kramer 
AMXSY-GA,  W.  Brooks 
AMXSY-J,  A.  LaG range 
AMXSY-L,  J.  McCarthy 
AMXSY-RA,  R.  Scungio 

M.  Smith 

3  Cdr,  USATECOM 
ATTN:  AMSTE-SI-F 

AMSTE-LFT,  D.  Gross 

R.  Harrington 

1  Cdr,  USACSTA 
ATTN:  STECS 


Intentionally  left  blank. 
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I 


USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes. 
Your  comments/answers  to  the  items/questions  below  will  aid  us  in  our  efforts. 

1.  BRL  Report  Number _ BRL-TR-3130 _ Date  of  Report  AUGUST  1990 

2.  Date  Report  Received _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest 

for  which  the  report  will  be  used.) _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source 
of  ideas,  etc.) _ _ 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars 
saved,  operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 

- $ - 

& 

.  satjszstsef  “  “ 

■  •  *6.*  'GeaasJ  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate 

-  changes?®  organization,  technical  content,  format,  etc.)  . . 

-  _ 

-  -  *  ’  - - - “ - 

*  -a  .TCaassaa - - - - — - 

Name 


CURRENT  Organization 

ADDRESS  _ 

Address 


City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  New  or  Correct 
Address  in  Block  6  above  and  the  Old  or  Incorrect  address  below. 


Name 


OLD  Organization 

ADDRESS  _ 

Address 


City,  State,  Zip  Code 


(Remove  this  sheet,  fold  as  indicated,  staple  or  tape  closed,  and  mail.) 


- FOLD  HERE - 


Department  of  the  Army 

Director 

U.S.  Army  Ballistic  Research  Laboratory 

ATTN:  SLCBR-DD-T 

Aberdeen  Proving  Ground,  MD  21  O  ' -5066 

OFFICIAL  BUSINESS 


BUSINESS  REPLY  MAIL 

FIRST  CLASS  PERMIT  No  0001,  APG,  MD 


POSTAGE  WILL  BE  PAID  BY  ADDRESSEE 


Director 

U.S.  Army  Ballistic  Research  Laboratory 

ATTN:  SLCBR-DD-T 

Aberdeen  Proving  Ground,  MD  21005-9989 


FOLD  HERE 


